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The  objective  of  this  report  is  to  theoretically  investigate  the  potential 
of  several  Igw  profile  antennas  for  use  at  or  near  ground  level  locations,  as 
compared  with  a  quarter  wavelength  monopole.  Of  specific  interest  is  the  radiated 
field  intensity  at  low  elevation  angles  produced  by  VHF  multitum  loop  antennas 
having  various  orientations  with  respect  to  a  flat  lossy  earth. 

To  accoiiplish  the  objective  stated  above,  we  start  with  Norton's  classical 
expressions  for  an  antenna  over  a  flat  lossy  earth  and  obtain  a  qualitative  indi¬ 
cation  as  to  how  the  presence  of  a  lossy  earth  affects  the  far- field  power  pattern 
of  an  antenna  in  air  located  above  the  lossy  °arth.  This  includes  a  description  of 
the  fields  in  terms  of  space  waves  and  surface  waves,  and  a  parameter  study  to 
determine  the  effects  on  the  far-field  pattern  of  varying  the  pattern  radius,  the 
ground  constants,  and  the  antenna  height  above  the  earth.  Next  a  conputati onal 1y 
straightforward  method  for  determining  the  far-field  pattern  of  an  antenna  located 
near  a  lossy  earth  is  described.  This  method  is  then  used  to  calculate  the  far- 
field  power  patterns  of  a  few  electrically  small  multitum  loop  antennas. 
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The  objective  of  this  report  is  to  theoretically  investigate 
the  potential  of  several  low  profile  antennas  for  use  at  or  near 
ground  level  locations,  as  conpared  with  a  quarter  wavelength  monopole. 

Of  specific  interest  is  the  radiated  field  intensity  at  low  elevation 
angles  produced,  by  VHP  muiti turn  loop  antennas  having  various  orientations 
with  respect  to  a  flat  lossy  earth. 

To  acconpl isn  the  objective  stated  above,  we  start  with  Norton's 
classical  expressions  for  an  antenna  over  a  flat  lossy  earth  ar.d  obtain 
a  qualitative  indication  as  to  how  the  presence  of  a  lossy  earth  affects 
the  far- field  power  pattern  of  an  antenna  in  air  located  above  the  lossy 
earth.  This  includes  a  description  of  the  fields  in  terms  of  soace  waves 
and  surface  waves,  and  a  parameter  study  to  determine  the  effects  on 
the  far- field  pattern  of  varying  the  pattern  radius,  the  ground 
constants,  and  the  antenna  height  above  the  earth.  Next  a  compu¬ 
tationally  straightforward  method  for  determining  the  far-field 
pattern  of  an  antenna  located  near  a  lossy  earth  is  described.  This 
method  is  then  used  to  calculate  the  far-field  power  patterns  of  a 
few  electrically  snail  multi  turn  loop  antennas. 


AN  INVESTIGATION  Of  THE  EFFECTS  OF  A  LC: 
EARTH  ON  ANTENNA  PATTERNS  AT  VHF 


I.  INTRODUCTION 

The  objective  of  this  report  is  to  theoretically  investigate 
the  potential  of  several  low  profile  antennas  for  use  at  or  near 
ground  level  locations  such  that  the  low  angle  radiation  can  be  com¬ 
pared  with,  that  of  a  quarter  wavelength  monopole  naving  tout  quarter 
■wavelengths  radial  arms.  Since  the  monopole  is  not  necessarily  vertical 
with  respect  to  the  earth,  we  are  interested  in  pattern  results  for 
various  tilt  angles  of  both  the  monopole  and  the  low  profile  antennas 
under  investigation. 

The  principle  low  profile  antenna  of  interest  is  the  multi  turn 
loop  antenna  on  a  small  ground  plane  tilted  at  various  angles  with 
respect  to  the  lossy  earth.  Thus,  there  is  a  need  to  develops  a  method 
for  calculating  the  far-field  pattern  of  an  arbitrarily  shaped  antenna 
having  an  arbitrary  orientation  with  respect  to  the  earth.  Therefore, 
a  secondary  objective  of  this  report  is  t.o  describe  a  method  for  doing 
this. 


In  Section  II  a  computationally  straightforward  method  for 
performing  the  calculations  described  above  is  presented.  Tne  method 
is  based  on  the  assumption  that  the  antenna  current  is  not  affected 
by  the  lossy  ground  and  gives  accurate  radiation  patterns  for  certain 
antennas  over  a  flat  lossy  earth  void  of  protruding  vegetation, 
buildings,  etc.  In  Section  II,  then,  are  presented  the  field  equations 
used  in  determining  the  antenna  patterns.  The  equations  are  interpreted 
in  terms  of  space  waves  and  surface  waves 

In  Section  III  a  parameter  study  is  made  using  the  theoretical 
expression  of  Section  II  to  determine  qualitatively  the  effects  on  an 
antenna's  far-field  pattern  of  the  distance  from  the  antenna  to  the 
observation  point,  the  antenna  geometry,  and  the  conductivity  and 
permittivity  of  the  ground. 

In  Section  IV  we  show  how  a  continuous  current  distribution  can 
be  expressed  in  terms  of  a  finite  number  of  infinitesimal  sources, 
and  in  Section  V  the  results  of  Sections  II  and  IV  are  used  to 
calculate  the  far-field  power  patterns  of  four  antennas. 

Appendix  I  contains  descriptions  and  listings  of  computer 
programs  used,  and  Appendix  II  compares  patterns  calculated  by  the 
methods  of  this  report  with  those  found  elsewhere  in  the  literature. 
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II.  iiilOi'Y 

In  this,  st'ction  wr  wish  to  drscriln*  ,j  niiijpijS.it ion.ii  ly  slr.ii«|hl 
forward  method  for  detemrininy  the  radiated  lit'lds  of  an  antenna 
located  over  an  infinite  flat  ground  plane  of  arbitrary  permittivity, 
and  conductivity.  A  number  of  approximations  are  made  and  these  are 
described  along  with  their  resulting  limitations  in  the  text  to 
follow.  The  method  is  best  applicable  to  electrically  small  antennas 
which  are  at  least  partially  shielded  from  the  lossy  earth  by  a 
small  perfectly  conducting  ground  plane  or  by  a  cavity,  and  yields 
the  radiated  fields  to  witnin  a  few  wavelengths  of  the  antenna. 
Electrically  large  antennas  can  be  treated  through  the  use  of  increased 
computer  time. 

The  basic  approximation  in  this  method  is  that  the  antenna 
current  with  the  antenna  located  in  free  space  be  identical  to  the 
current  with  the  antenna  located  over  a  finitely  conducting  flat 
earth.  Rigorously  this  is  clearly  not  the  case.  However,  it  is  a 
reasonable  approximation  for  an  antenna  which  is  shielded  from 
the  finitely  conducting  earth  as  is  the  case  for  an  antenna  in  a  - 
cavity  or  over  a  small  ground  plane,  and  it  is  these  configurations 
that  we  will  consider.  A  typical  configuration  is  Shown  in  fig.  1. 
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FINITELY  CONDUCTING  GROUND 


rig.  1.  Antenna  shielded  from  ground  by  perfectly 
conducting  ground  plane. 


Since  small  changes  in  a  current  distribution  cause  even  smaller 
changes  in  the  far- field  pattern,  our  approximation  is  most  applicable 
to  determining  far-field  Dattems. 

The  problem  of  determining  the  far-field  pattern  of  an  antenna 
must  begin  with  a  determination  of  the  currents  on  the  antenna.  When 
the  antenna  is  in  the  presence  of  a  conducting  earth  this  can  be 
done  using  the  rigorous  Sonmerfeld-integral  approach[l-3],  or  by 
some  approximate  technique,  one  of  which  is  a  reflection  coefficient 
approximation  described  by  E.  K.  Miller  et  al[4,5].  The  reflection 
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coefficient  approximation  is  a  technique  using  the  Fresnel  plane- 
wave  reflection  coefficients  to  account  for  the  influence  of  fields 
reflected  from  the  interface  upon  the  antenna  current  distribution, 
and  is  found  to  yield  the  input  impedance  generally  to  within  10‘* 
of  the  results  obtained  by  the  rigorous  method  but  at  about  12  of 
the  computer  cost.  Our  approach  is  more  approximate  than  either 
of  the  above,  but  has  the  advantage  that  computer  programs  already 
exist  to  evaluate  the  currents  on  antennas  located  in  free  space[6-9j. 
As  previously  mentioned,  it  is  most  applicable  when  the  antenna  is 
somehow  shielded  from  the  earth.  However,  it  would  also  apply  to 
antennas  sufficiently  elevated  from  the  earth. 


Assume  now  that  by  some  method  a  current  distribution  has  been 
determined  on  an  antenna  or  a  wire-grid  model  of  the  antenna.  To  find 
fhc-  fields  radiated  by  this  arbitrary  current  distribution  we  must 
first  find  the  fields  radiated  by  an  arbitrarily  orientated  infinitesimal 
electric  current  element  located  in  air  and  near  a  plane  earth  of 
relative  permittivity  er,  conductivity  o  expresses  in  electromagnetic 
units,  and  permeability  u0  of  free  space.  (Note  one  electromagnetic  unit 
or  emu  equals  10"*'  mho/neter).  To  do  mis  we  can  use  expressions  given 
by  nortonpO], 

Figure  2  shows  the  coordinate  system  to  be  used  We  let  the 


z 


i,  j,  and  k  in  the  direction  of  x,  y,  and  z.  The  cylindrical 
coordinates  r,  c,  and  z,  and  corresponding  unit  vectors  r,  o,  and  k 
are  also  shown.  Ri ,  R2,  and  a  are  defined  in  Fig.  2  and  in  addition 

(1 )  k  2t,/i  ,  a  =  free  space  wavelength 

(2)  =  k2(tr-jx) 

l^j\  _  T  f:  ..  1/-.18  /r 


=  1.8  X  10'°  3  /f, 

emu'  kc 

-  ;.ojnd  conductivity  eAp  ssed  in  emu 


f,  =  frequency  in  kilocycles  per  second 

a2  =  x2+y2+z2 


=  k/k. 


r  Norton  found  simplified  and  computationally  useful  expressions 
for  E  and  we  will  include  some  of  them  here.  While  Norton's  expressions 
are  valid  in  the  upper  half-space  to  within  a  few  wavelengths  of  the 
source,  the  field  components  to  follow  are  simplified  for  the  far-field. 
By  far-field  in  this  report  we  mean  1/R  »  l/R^.  For  a  vertical  electric 
current  element,  with  time  dependence  surpressed,  and  located  at 

z=a,  x=y=0,  as  shown  in  Fig.  2, 


-  ik  -< 


jk  -jcosV 


+  R  COS2V*  “r 

v  ■- 


-jkR2  1 

+  (1-kv)(1-u2+u4  COS2--;' )  F  ^ - 

2  i 


where. 

(9a) 

R» 

=  (sino'  -  ujl-u2cos2*'  )/(sin-j 

(9b) 

sin?' 

=  (z+a)/R2 

(9c) 

F 

=  [1-j  fnwe~w  erfc(j*/w)j 

(9d) 

w 

=  4Pl/(l-Rv)2 

(9e) 

h 

=  p  e‘Jk 

rR0  cos  b 

(9f) 

p 
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(Sg) 

(9h) 


tan  b  =  Ur  +  cos2y'  }/x 

sin/'  =  (z-a)/R-|  . 

In  Eq.  (8)  the  first  two  terms  represent  the  direct  or  space  wave 
radiation  from  the  source  and  its  image  (located  at  z=-a,  x=y=0) 
respectively,  and  both  vary  as  I/R.  The  last  term  represents  the 
ground  or  surface  wave  radiation  and  varies  as  F/R.  F  is  generally 
termed  a  ground  wave  attenuation  function,  and  depends  upon  the 
earth's  constants  and  upon  the  distance  to  the  field  point,  p  and  b 
are  generally  termed  the  numerical  distance  and  phase  factor 
respectively  for  vertical  dipoles.  The  remaining  components  for  the 
vertical  element  are 
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-jkRn 


=  jk  -;sino“cos<"  e 


+  R  sin/cos.;. 
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-  coso'  (1 -Ry  )u 'll -u^ cosin'  F  ~ 


-jkR, 


2 


L  u2(l-u2cos2y' )  sin2-/  \  : 

v  2  *  2  )  i 


jkRp 


"2 


(ID 


Ev  0. 


For  a  horizontal  infinitesimal  electric  current  element  aligned 
with  the  x-axis 
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-jkR 


=  -  jkcosMSin/'coso 


. »  e 
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-jkR, 


-  K  siny  cos=, 
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+  COS; 


(1  -R  )u  jl-U2COS2*'  F  ^ 


-jkR- 


'1 


03) 


f,  u2d-u2cosV)  si  ^  I 

'  V  2  *-J{ 


:-‘r  =  jkcoso  s’ sin2/'  ^ 


-jkR 
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1 

2  e  c 
R  si nV  - - 

V  K„ 


+  d+Rh)  G 


-jkR, 

- cdsW{1-R  )  F 

'2  v 


-jkR,  } 
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(14) 

where. 


Eh  = 


-  jksin<j 


:  -jkR-j  -jkR2 


-  jk 


-  R, 


i 

'v 


1 


h  R„ 


+  O+Rh) 


(15a) 
(15b) 
(15c) 
(15d) 
(1 5e ) 
(15f) 


R.  =  ( i1-u2cos2i‘  -  usin-v*  )/(f1-u2cos2v'  +  usin-/) 
n 

6  =  [1-j/nv  e“v  erfc(j/v)] 
v  =  4qi/(HRh)2 
q,  =  -q  ejV 

«  n 

*i,  "  'f'r'2 

^  "  >.cos  (b* ) 

tan(h’)  =  (e^-cos2*’ )/x  * 


G  is  a  ground  wave  attenuation  function  similar  to  F.  q  and  b‘  are 
generally  termed  the  numerical  distance  and  phase  factor  respectively 
for  horizontal  dipoles.  In  Eqs.  (8)  and  (10 )- (14 )  surface  waves 
terms  can  be  identified  by  the  presence  of  F  or  G  as  a  facto.'.  Eqs. 
(8)  and  nO)-(ll)  including  higher  order  terms  were  programmed  for 
the  Datacraft  6024  digital  computer.  Subroutine  DIPOLE  is  described 
i n  Appendi x  I . 

We  now  wish  to  investigate  the  nature  of  the  ground  wave 
attenuation  functions.  Since  they  are  similar  we  will  look  only 
at  F.  Combining  Eqs.  (9d-f)  we  have 

*2 


If  we  consider  a  point  in  the  far-field  and  at  a  low  elevation  angle, 
such  that  R2  •»  (z+a),  then  from  Eqs.  (9a,b)  we  have  Ry  %  -1  and 

R 

(17)  w  %  -cos(b)  e~^b  • 

Inserting  the  value  of  *  from  Eq.  (3)  into  Eq.  (17)  we  have 

p;  f 

(18)  W  £  -cos(b)  e”jb  — - — 

1 .8x10  ao 


(16) 


w 


4ncos (b)  e  J 


ib 


(  <- 


k  r 

v 
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For  frequencies  near  200  \iz,  and  a  typical  ground  conductivity 
of  10“'3  emu,  w  will  be  on  t'ne  o  Jer  of  R2  (meters)  in  magnitude. 

From  Eq.  (9g)  we  see  0  ^  b  ^  r,/Z  and  thus  — /2  arg  w  _  0.  Thus, 
the  argument  of  the  complementary  error  function  of  Eq.  (9c)  is  a 
complex  nimber  whose  phase  is  between  r,/ 4  and  -r/2 .  Its  magnitude  may 
be  evaluated  exactly  using  Eq.  (13),  but  for  frequencies  around 
200  MHz  and  for  typical  ground  constants  (i.e.,  e*.  %  5,  j  %  10~^3  emu) 
its  magnitude  is  on  the  order  of  ^  £  &  (meters/.  In  this  case 
the  following  asymptotic  expansion  for  erfc(z  =  jVT)  is  validfllj 


09) 


,2 

»T  z  e  erfc(z)  -v  1  +  ) 

m=l 


1-3 —  {».- 
(2 z2)” 


I  \ 


Retaining  only  the  first  term  of  the  summation  we  have 

72  9 

(20)  /T  z  ez  erfc(z)  ^  1  +  l/2z  . 


If  in  Eq.  (9c)  we  make  the  change  cf  variables  =  z  we  get 

_,2 

(21)  F  =  [1  -  z  e  z  erfc(z)]. 

Combining  Eqs.  (20)  and  (21)  v/e  get  finally  the  well  known  result 

(22)  F  -v  -  l/2z^  %  constant/R. 


Since  surface  wave  terms  vary  as  F/R  or  G/R  we  then  see  that  they  have 
a  1/R2  distance  dependence  for  VHF  frequencies  and  typical  ground 
constants. 

We  next,  wish  to  consider  tne  far- field  space  wave  at  low 
elevation  angles.  Noting  that  cosv"  %  cosv’  %  1,  and  P.  %  -1,  the 
space  wave  portion  of  Eg.  (8)  becomes  v 

/-jkP^  -jkR^ 

(23)  E^(space)  =  -  jk  - ^75 - -  j  . 

z  ^  "2  J 

We  can  write  in  terms  of  R-j  using  the  approxi -nation 

(24)  +  aR  =  +  2asina  %  R^+2aa. 

Combining  Eqs.  (23)  and  (24)  and  using  1/R^  R;  l/R^  %  1  /R  we  get 
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(25) 


__  j  ^ 

f-2  (space)  %  2ak  e  1  (|)  . 


Thus,  for  a  given  elevation  angle  u  the  space  wave  has  the  usual 
1/R  dependence,  and  we  would  expect  it  to  predominate  the  surface 
wave  except  at  elevation  angles  so  small  that  a  is  on  the  order  of 
or  smaller  than  1/R.  Eq.  (25)  also  predicts  that  the  space  wave  field 
strength  should  increase  with  increasing  elevation  angle  or  increased 
height  of  the  vertical  dipcle  from  the  earth's  surface.  Might  on  the 
earth's  surface,  where  -  0,  the-  space  wave  vanishes  and  ail  radiation 
is  due  to  the  surface  wave.  Summarizing  we  see  that  the  surface  wave 
predominates  at  very  low  elevation  angles,  and  the  space  wave  at  the 
higher  elevation  angles. 


From  the  above  discussion  we  do  not  mean  to  imply  that  the 
peak  magnitude  of  the  surface  wave  field  is  always  considerably  less 
than  the  peak  magnitude  of  the  space  wave  field.  The  numerical  distance 
can  be  related  to  or  -Gj  and  can  thus  be  a  measure  of  the 
relative  importance  of  the  surface  wave.  In  i'.g.  ;  is  plotted  Mif* 


Fig.  3.  The  absolute  value  of  the  ground  wave  attenuation 
function  for  vertical  dipoles  versus  the  numerical 
distance. 


absolute  value  of  the  ground  wave  attenuation  function,  IF;,  versus 
p,  the  numerical  distance  for  vertical  dipoles,  and  we  see  that  for 
p  -  10,  ?Fi%  l/2p[14].  Thus  the  surface  wave  is  less  important  for 
larger  values  of  p.  If  we  consider  the  exarrple 

f  =  200  MHz, 

R  =  1  km,  and 
o  =  10"^  emu. 


then  we  have  p  £  250  cos (b) ,  and  from  Fig.  3  F;  -c  0.002.  However,  if 
we  consider  the  same  example  but  witn  f  =  I  KHz  we  have  p  ^  .005  cos(b} 
and  |F'  £  1.0.  The  horizontally  polarized  surface  wave  is  generally 
attenuated  more  rapidly  than  the  vertical  polarization.  Jordon  and 
Balmain[14]  present  a  good  discussion  of  the  effects  of  the  ground 
constants,  frequency,  and  range  on  the  pattern  of  infinitesimal  electric 
current  elements. 


III.  NUMERICAL  PARAMETER  STUDY 
A.  Introduction 


In  this  section  we  wish  to  look  at  the  far- field  power  oattems 
of  the  basic  radiating  elements,  which  are  the  horizontal  and  vertical 
infinitesimal  electric  current  elements.  The  element  patterns  were 
calculated  for  the  planes  shown  in  Fig.  4.  Due  tj  symmetry  only  the 
first  90°  of  the  pattern  in  the  xz  plane  will  be  shown  for  the  vertical 
element.  Patterns  for  the  horizontal  element,  assumed  to  be  aligned 
with  the  x  axis,  will  be  shown  for  the  first  90°  of  the  plane  containing 
the  dipole  axis  and  the  plane  perpendicular  to  the  dipole  axis,  i  e, , 
the  xz  and  y z  planes  respectively. 

The  pattern  for  an  antenna  over  a  lossy  earth  is  dependent  upon 
many  factors.  Sone  of  the  most  inportant  are: 

1.  Radial  distance  from  the  source. 

2.  Geometry  of  the  antenna,  and  the  antenna  height  from  the 
surface. 

5.  Conductivity  and  permittivity  of  the  lossy  ground. 

4.  In  a  practical  sense,  the  vegetation,  buildings,  vehicles, 
etc. 

Clearly  it  is  impractical  to  show  the  influence  of  the  first  three 
parameters  each  time  an  antenna  pattern  is  displayed  and  virtually 
impossible  to  account  for  the  fourth.  Thus,  we  wish  to  do  a  parameter 
study  on  the  infinitesimal  electric  current  element  to  obtain  a 
qualitative  indication  for  the  effects  of  the  first  three  parameters 
on  real  antennas.  Buildings,  vegetation,  hills,  and  other  obstacles 
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Fig.  4.  Geometry  for  patterns. 


effect  an  antenna's  pattern  in  two  ways.  They  may  either  scatter 
or  absorb  some  of  the  energy  radiated  by  the  antenna,  and  this 
scattering  or  absorption  may  have  a  much  greater  effect  on  an 
antenna's  pattern  than  any  or  all  of  the  first  three  parameters.  Not 
knowing  the  geometry  of  the  obstacles  one  can  only  assume  that  the 
amount  of  radiation  received  at  one  point  from  a  source  located  at 
another  point  is  directly  proportional  to  the  amount  of  radiation 
received  in  the  absence  of  the  obstacles.  This  is  our  justification 
for  studying  the  problem  of  propagation  over  an  infinite  flat  lossy 
earth.  Stated  another  way,  we  wish  to  show  the  effect  of  the  lossy 
earth  on  an  antenna's  pattern,  while  we  understand  that  in  practice 
there  are  other  factors,  (i.e.,  scattering  or  absorbing  obstacles) 
affecting  the  pattern.  Unless  it  is  one  of  the  parameters  being 
varied,  patterns  will  be  shown  for  R  -  1  km,  cr  =  5,  t  =  3  x  10"'3 
emu  =  0.03  niio/meter,  and  f  =  162  MHz.  A  log  log  scale  was  chosen  to 
display  this  data  because  it  enphasizes  the  low  elevation  angles  of 
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interest  in  this  study  and  provides  for  a  good  dynamic  range.  For  each 
of  tiie  three  parameter  studies  the  patterns  are  normalized  seen  that  the 
largest  peak  power  for  any  curve  shown,  within  a  particular  parameter 
study,  is  unity.  Thus,  even  though  a  particular  study  may  contain  marry 
graphs,  each  having  many  curves,  any  two  curves  may  be  compared  directly. 

The  pattern  calculations  in  the  following  sections  were  made 
using  subroutine  DIPOLE  described  in  Appendix  I.  Subroutine  DIPOLE 
includes  the  effects  of  space  waves  and  surface  waves,  even  tnougn  fer 
many  of  the  calculations  only  the  space  wave  is  significant.  In  these 
cases  much  effort  may  be  saved  by  using  only  the  first  two  terms  in 
Eqs.  (8)  and  (10;-(14). 

B.  Pattern  Radius 

In  this  and  in  following  sections  various  far-field  power  patterns 
will  be  given  for  R  =  1  km.  If  these  patterns  are  assumed  to  apply 
to  other  radii  by  simply  introducing  the  usual  1/3^  power  dependence 
some  error  will  be  incurred,  especially  for  low  elevation  angles. 

Using  the  1 /r2  dependence  essentially  amounts  to  ignoring  the  surface 
wave  and  assuming  the  space  wave  predominates.  As  was  seen  in 
Section  II,  this  approximation  will  be  valid  except  at  low  elevation 
angles.  Figures  5,  6  and  7  have  been  drawn  to  obtain  an  idea  of  the 
nature  of  this  error.  Ihe  two  curves  of  Fig.  5  are  the  relative  errors 
in  extrapolating  from  1  km  to  0.1  km  or  1  km  to  10  km  using  a  1/R2 
power  dependence  for  a  vertical  dipole  located  */4  above  a  ground 
with  cr  *  5,  o  =  .03  mho/meter  -  3  x  10-^3  emu,  and  f  =  162  KHz.  If 
Pa  is  the  power  at  Ra  and  some  elevation  angle,  and  Pfa  is  the  Dower 
at  R5  and  the  same  elevation  angle,  then  the  relative  error  in 
extrapolating  from.  Pa  to  f%  is  given  by 


Figures  5-7  were  drawn  by  calculating  the  error  exactly  (i.e., 
using  subroutine  CIRCLE)  every  0.1°  from  0.1°  to  1.0s  and  every  degree 
from  1®  to  10°.  The  error  at  any  point  is  not  only  a  function  of  the 
relative  magnitude  of  tie  space  and  surface  waves,  but  also  of  their 
relative  phases.  Thus,  .it  a  point  where  the  space  and  surface  waves 
add  in  the  proper  phase  the  error  ray  be  considerably  less  than  at 
nearby  points.  This  accounts  for  the  sometimes  erratic  nature  of  the 
curves  in  Figs.  5-7. 

In  Fig.  5  we  see  that  it  is  possible  to  extrapolate  from  1  to 
0.1  km  with  a  relative  error  less  than  IS  for  elevation  angles  greater 
than  6°,  while  extrapolation  from  1  to  10  km  can  be  done  with  a  relative 
error  less  than  1*  for  elevation  angles  greater  than  2°.  Figures  6  and 
7  are  identical  to  Fig.  5  except  that  they  are  for  horizontal  dipoles 
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with  the  pattern  in  the  xz  and  yz  planes  respectively.  Note  in 
Figs.  5-7  that  extrapolation  to  greater  radii  is  done  more  accurately 
than  to  shorter  radii.  This  is  because  the  surface  wave,  which  is 
the  cause  of  the  error,  is  less  significant  at  greater  distances. 

C.  Dipole  Orientation  and  Height  from  Surface 

Here  we  wish  to  obtain  an  indication  of  the  type  of  pattern 
radiated  by  the  infinitesimal  vertical  dipole  and  the  horizontal 
dipole  in  either  plane.  We  also  desire  to  determine  the  effect  of 
elevating  the  dipoles  from  tr.e  earth's  surface. 

In  Fig.  8  curves  are  shown  for  the  vertical  element,  and  Figs. 

9  and  10  contain  curves  for  the  horizontal  element  in  the  xz  and  yz 
planes  respectively.  We  note  the  following  from  Figs.  8-10. 

1.  As  one  would  expect,  the  vertical  elements  are  the  best 
radiators  at  low  elevation  angles,  and  this  is  especially  true  when 
the  elements  are  less  than  a/2  from  the  earth's  surface.  The 
horizontal  elements  radiate  best  at  very  large  elevation  angles. 

2.  Except  for  the  horizontal  dipole  in  the  xz  plane,  which  is 
the  poorest  radiator  at  low  elevation  angles,  raising  the  element 
from  the  earth’s  surface  greatly  increases  the  amount  of  radiation 
at  low  elevation  angles. 

3.  The  radiation  from  an  elevated  element  is  the  sum  of 
radiation  from  the  source  and  its  image.  Thus,  raising  the  element 
makes  the  source  look  larger  and  sidelobes  appear  in  the  patterns. 

In  Fig.  11  we  plot  the  elevation  angle  corresponding  to  the 
first  pattern  maximum  versus  the  dipole  height.  Note  that  for  the 
vertical  element,  and  the  horizontal  element  in  the  yz  pla-e  it  is 
possible  to  shift  this  maximum  well  under  2'’  by  elevating  8\  or  more, 
wnile  the  peak  of  the  pattern  for  the  horizontal  element  in  the 
xz  plane  is  never  less  than  40°. 

Figures  8-10  can  be  used  as  qualitative  design  information 
since  the  fields  of  an  extended  source  is  simply  an  integration  or’ 
sumration  of  the  fields  of  infinitesimal  sources.  For  example,  if 
one  wants  to  design  an  antenna  to  radiate  at  low  elevation  angles, 
and  is  restricted  to  placing  his  antenna  very  near  the  earth's  surface, 
then  he  should  rely  on  vertical  wires.  However,  if  radiation  at  high 
elevation  angles  is  desired,  then  horizontal  wires  should  be  used. 

As  a  second  example  assume  one  has  the  pattern  for  a  horizontal 
linear  antenna  located  A/2  above  the  earth  in  the  plane  containing 
the  antenna  axis.  It  is  desired  to  find  the  effect  of  lowering  the 
antenna  almost  to  the  earth's  surface.  From  Fig.  10  we  see  that  this 
■will  seriously  degrace  the  antenna's  performance,  especially  at  the 
lower  elevation  angles. 
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Fig.  11.  Angle  of  pattern  maximum  nearest  the  ground 
versus  the  dipole  height. 


L1.  Ground  Constants 

Here  we  wish  to  determine  the  effects  on  a  power  pattern  cf 
varying  the  permi tti vi ty  and  the  conductivity  of  the  earth.  The  range 
of  ground  constants  covered  is  .  r  =  2  to  32  and  ■  =  lO-'^  to  10"^J 
emu  or  0.1  to  0.0001  mho/meter.  All  curves  are  for  the  infinitesimal 
dipoles  located  x/4  above  the  earth's  surface  and  at  P.  =  1  km. 

Figure  12  is  for  the  specific  case  or  a  vertical  dipole,  ?  =  1 0-^ 2  emU} 
and  curves  are  drawn  for  zr  =  2,  4,  8,  16,  and  32.  Figures  1 3-15 
are  identical  to  Fig.  12  except  that  a  =  10"^3}  lo-1^,  and  1G-^5  emu 
respectively.  Figures  16-19  are  identical  to  Figs.  12-15  respectively 
except  that  they  are  for  a  horizontal  dipole  with  the  pattern  in 
the  xz  plane,  and  Figs.  20-23  are  identical  to  Figs.  16-19  except 
that  they  are  fcr  a  horizontal  dipole  in  the  yz  plane. 

Note  that  in  Figs.  12-23  varying  c  over  3  orders  of  magnitude 
has  almost  no  effect  on  the  patterns,  while  varying  er  from  2  to  32 
does  significantly  alter  the  patterns  especially  for  the  vertical 
dipole.  Increasing  cr  increases  the  radiation  at  low  elevation 
angles  for  vertical  dipoles,  and  decreases  the  radiation  at  low 
elevation  for  horizontal  dipoles.  The  radiation  from  horizontal 
dipoles  at  large  elevation  angles  is  increased  as  > r  is  increased. 

Figures  12-23  could  be  used  if  one  knew  for  example  the  pattern 
of  a  vertical  linear  antenna  over  a  grourd  with  %  4  and  o  -l  10'^ 
emu,  and  it  was  desired  to  qualitatively  predict  its  performance  if  it 
were  located  ever  a  ground  with  er  %  10  and  %  IQ-'5  emu.  Then 
using  Figs.  12  and  14  we  see  that  the  power  at  low  elevation  angles 
will  be  increased  slightly,  while  that  at  large  elevation  angles  will 
be  essentially  unchanged. 


IV.  REPRESENTING  COfiTI TiLiCLS  CURRENT  DISTRIBUTIONS 
BY  A  FINITE  NUMBER  OF  INFINITESIMAL  SOURCES 


In  this  section  we  wish  to  devise  a  method  for  determining  the 
fields  radiated  by  an  arbitrary  current  distribution  in  terms  cf  & 
finite  summation  cf  the  fields  of  infinitesimal  sources  rather  than 
by  an  integration  over  the  source  distribution.  Consider  the  wire 
in  Fig.  24a  which  is  assui.°d  to  carry  the  arbitrary  current  distribution 
T($').  The  fields  of  this  current  distribution  at  the  field  point  f; 
can  be  expressed  using  the  following  linear  vector  function 


(27) 


? 


wire 


20 


Patterns  for  vertical  dipole  with 


"  vertical  dipole  with  a=10  emu  versus 


ELEVATION  ANGLE  (DEGREES) 
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ELEVATION  ANCLE  (DEGREES) 
Patterns  i"or  horizontal  dipole 


ION  ANGLE  (DEGREES) 


where. 


(28)  A  ;  ir  E,ncosc  +  i?E,Jlsinc  +  ikE'^cos^ 
r  ■?  2 

+  jrE'^sin?  -  jcE'^cos-?  +  jkE'^sin? 
r  o  4 


+  krE' 
r 


+  kkE^. 


Note  that  the  jo  term  of  A  differs  in  sign  from  that  given  by  NcrtonpO], 
and  it  is  felt  that  this  is  probably  a  typographical  error  ir.  Norton's 
paper. 

The  unprimed  electric  field  components  are  defined  in  Eos.  (8) 
and  (10)-(14),  and  the  primed  components  are  related  to  the  unprimed 
components  by 


(29a) 

(29b) 

(29c) 


-h 


r-.h 

c  rcos$ 


E*;  =  E'hsin6 

6  9 

=  E'^cos? 


o  and  «i  are  defined  in  Fig.  2  and  are  in  general  functions  of  (R-R * ) . 
Comparing  Eqs.  (10)- (14 )  and  (28),  (29)  we  see  that  ail  field  comoonents 
in  Eq.  (28)  can  be  evaluated  without'  dividing  by  sine  or  cose. 

It  seems  reasonable  that  the  fields  radiated  by  I(R')  would 
be  closely  approximated  by  the  fields  radiated  by  several  inf'  tesimal 
current  elements  orientated  col  linear  with  and  located  along  uie  wire, 
and  having  a  complex  magnitude  related  to  the  complex  magnitude  of 
I(R')  at  their  locations.  To  make  this  more  clear  consider  for 
example  a  vertical  linear  antenna  of  length  L  located  along  the  z  axis 
as  shown  in  Fig.  24b.  In  this  case  Eq.  (27)  becomes 


(30) 


E(R)  = 


■> 

-*■ 


I (z 1 )  *  A(R-z')dz'. 


Our  approach  will  be  to  numerically  evaluate  the  integral  of  Eq.  (30). 
For  example  we  could  use  the  trapezoidal  rule  which  is  the  approxi¬ 
mation 
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(31) 


X  —  X 

f (x)dx%  2y -[ t'Cx-j  hf(x2)j. 


In  this  case  Eq.  (30)  becomes 
(32) 
where 


»  ->  !  -* 

»■  /  ft  \  r,  ^  r  f  /  _  i  '  .  ^  i  ' 

<  2^-i- \Za;~rv z^/ j 


F(2' )=I (2* )*A(R-2* ) . 


In  essence  we  are  approximating  the  continuous  current  distribution 
of  Fig.  23b  by  two  infinitesimal  current  elements  located  at  z‘a  and  z£ 

j  -*» 

and  of  magnitude  and  orientation  given  by  jl(za)  and  (2^ }  respectively. 

lit'  there  is  considerable  variation  of  the  current  over  the  length  of 
tne  wire  we  may  wish  to  sol it  the  integral  of  Eq.  (30)  into  two  parts 
to  get 


(33) 


■  \  ,r/  a  b) 

’  "C  'V*F(-T-j 


% 


,  i-j.v  -*iK+zL> 

ijlr /z*  i+F  !  _2 _ tL\ 

Z\?  lV  V  2  j 


Clearly  we  could  continue  this  splitting  to  obtain 


(34) 


AL[lf(z;)+F(z|h---+F(z‘. 


where. 


(35a)  aL  =  (z'-z’)/N 

(35b)  2n  =  za  +  |  L  n=0,l ,2,- **N. 


Inspection  of  Eq.  (34)  shows  that  we  are  representing  the  continuous 
current  distribution  by  N+l  infinitesimal  current  elements  located 
at  Zp  and  with  complex  magnitude  and  orientation 
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otherwise. 
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Any  numerical  integration,  such  as  Gaussian  quadrature  or 
Simpson’s  rule,  will  permit  representing  a  continuous  distribution 
by  a  series  of  infinitesimal  current  elements  whose  location, 
orientation,  and  magnitude  are  most  obvious  after  the  numerical 
integration  is  written  out  explicitly. 

While  in  this  example  we  chose  our  current  distribution  to 
be  aligned  with  the  z-axis,  this  in  no  way  restricts  our  procedure. 
For  example  if  we  wish  to  approximate  the  fields  radiated  by  the 
current  distribution  in  Fig.  25  by  one  application  of  the  trapezoidal 
rule,  we  would  have 


i 
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Fig.  25.  Current  distribution  on  an  arbitrarily 
orientated  straight  wire. 


(36)  EXR)  =  J  I(R’)  -  A(R-R')  d&  ,Z;)*F(x£  ,y£ ,z’ )] 


wi  re 


where 


A  large  class  of  antennas  can  be  modeled  by  a  wire-grid  con¬ 
sisting  of  many  straight  wire  segments.  Often  several  segments 
connect  at  er,e  point.  In  representing  such  a  complex  structure  by 
a  finite  number  of  infinitesimal  electric  current  elements,  the 
path  of  integration  in  Eq.  (27}  is  over  the  entire  wire-grid 
structure.  To  avoid  errors,  it  is  recommended  that  this  integral  be 
split  into  smaller  parts  to  avoid  integrating  past  a  point  where  two 
or  more  wires  connect. 


Far- field  pa 
are  compared  with 


items  calculated  using  the  methods 
patterns  found  in  the  literature  in 


of  tuis  report 
Appendix  II. 


Thus  far  we  have  treated  only  current  distributions  $s  would 
exist  or.  a  wire  antenna,  or  if  a  wire-grid  model  was  made  of  a  complex 
antenna.  The  technique  can  however  be  extended  to  two  or  three 
dimensional  sources.  The  approach  would  be  to  numerically  evaluate 
the  surface  or  volume  integral  corresponding  to  Eq.  (k7).  This  would 
replace  the  surface  or  volume  integral  by  a  double  or  a  triple 
summation,  the  inspection  of  which  would  indicate  the  location  and 
magnitude  of  the  infinitesimal  electric  current  elements. 


We  will  now  summarize  the  ideas  of  the  preceding  section*. 

The  problem  we  wish  to  solve  is  that  of  finding  the  pattern  of  an 
antenna  in  the  vicinity  of  a  flat  finitely  conducting  earth.  We 
use  the  approximation  that  the  current  on  the  antenna  is  unaffected 
by  the  presence  of  the  earth,  and  thus,  may  be  calculated  assuming 
the  antenna  to  be  in  free  space.  This  is  a  reasonable  approximation 
if  the  antenna  is  shielded  from  the  earth  by,  for  example,  a  small 
perfectly  conducting  ground  plane  or  is  sufficiently  elevated  aoove 
the  earth.  Next  the  continuous  currents  on  the  wires  of  the  antenna 
are  represented  by  a  number  of  infinitesimal  current  elements.  The 
fields  of  the  individual  elements  can  be  found  using  Eqs.  (8)  and 
(10)-(14),  and  then  summed  to  yield  the  antenna  pattern.  This  pro¬ 
cedure  will  be  applied  to  several  antennas  in  the  next  section. 


V.  THEORETICAL  ANTENNA  PATTERNS 

We  now  wish  to  apply  the  methods  of  Sections  II  and  IV  to 
determine  the  far  field  patterns  of  four  antennas  of  particular 
interest: 

1.  A  quarter  wave  monopole  over  four  quarter  wave  arms  as 
shown  in  Fig.  26. 

2.  A  3-1/2  turn  multitum  loop  (KTL)  as  shown  in  Fig.  27. 

3.  A  12  turn  "picket  fence"  KTL  as  shown  in  Fig.  28. 
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Fig  26.  Antenna  1  -  A  quarter  wave  monopole  over  four 
quarter  wave  arms. 
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4.  A  12  turn  "picket  fence"  KTL  as  shown  in  Fig.  28,  but 
with  a  four  inch  diameter  reflecting  plate  placed  0.4"  above 
the  KTL . 


Antennas  2,  3,  and  4  are  placed  over  a  perfectly  conducting 
ground  plane  k/Z  in  diameter.  The  four  horizontal  quarter  wave  arms 
of  antenna  1  approximate  a  ground  plane  a/2  in  diameter.  The  small 
perfectly  conducting  ground  plane  is  then  placed  */3  over  a  lossy 
earth  with  cr  =  5  and  c  =  3.0  x  10“^  emu  =  0.03  mho/meter,  as  shown 
in  Fig.  28s.  Antenna  power  patterns  will  be  shown  for  R  =  1  kq, 
f  =  162  MHz,  and  for  various  antenna  tilt  angles  as  defined  in 
Fig.  29b.  For  large  tilt  angles  the  \/2  diameter  ground  plane  inter¬ 
sects  the  lossy  earth.  This,  a  problem  since  subroutine  DIPOLE, 
used  in  making  the  calculations  in  this  section,  applies  only  to 
sources  located  above  the  lossy  earth.  However,  since  we  are  mainly 
interested  in  seeing  trends  as  the  antenna  is  tilted,  this  pioblem 
is  avoided  by  simply  ignoring  the  few  sources  which  are  located 
beneath  the  earth's  surface. 

As  has  been  pointed  out,  a  calculation  of  the  far-field  pattern 
of  any  antenna  begins  with  a  calculation  or  an  approximation  of  the 
currents  flowing  on  the  antenna.  Thus,  for  antennas  2-4,  we  begin 
by  making  a  wi re-grid  model  of  the  MU  antenna,  the  >,/Z  diameter 
ground  plane,  and  for  antenna  4  the  4  inch  reflecting  plate.  Antenna 
1  is  identical  to  its  wire-grid  model.  Next  the  currents  flowing  on 
the  antenna  are  calculated  using  a  moment  method  solution,  and  the 
assumption  that  the  wire-grid  model  is  located  in  free  space.  As 
indicated  in  Section  IV  the  currents  flowing  on  the  wire-grid  model 
are  replaced  by  a  number  of  infinitesimal  electric  current  elements. 
Finally,  the  fields  of  these  infinitesimal  sources  are  calculated 
using  Eos.  (8),  (10)-(14),  and  (27),  (28)  and  added  vectorially  to 
form  the  antenna  field  patterns.  Power  patterns  are  obtained  from 
tiie  field  patterns  bv  taking  the  magnitude  squared  of  the  field  and 
multiplying  by  the  antenna  efficiency  defined  as 

F.  -  P,  P  . 

ir,n\  _  in  loss  _  rad 

{37}  Lfficiency  =  — y -  =  ~p — 

‘  i  n  i  n 

where  P.-n,  P-joss *  Prad  are  P0wer  input,  ohmic  power  losses,  and 
the  radiated  power  respectively.  The  efficiencies  were  measured[12] 
using  a  method  enployed  by  Wneeler[13j.  Tne  efficiency  of  a  ffTL  is 
strongly  dependent  upon  the  volume  it  occupies,  on  its  geometry,  and 
on  the  conductor  resistance  and  therefore  the  diameter  of  the  wire 
being  used.  A  study  is  presently  being  conducted  at  the  Electro- 
Science  Laboratory  to  optimize  the  efficiency  of  HTL  antennas.  At 
this  point  typical  results  are  35%  for  antenna  2,  and  20'%  for  antennas 
3  or  4.  However,  these  values  are  not  to  be  taken  as  being  the  best 
possible  for  a  g.ven  MTL,  out  sinply  typical  results  to  date.  For 
conparison,  the  efficiency  of  antenna  1  is  taken  to  be  70%. 
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LOSSY  GROUND 


(b) 


Fig.  29.  (a)  Geometry  cf  antennas  with  respect  to  the  lossy 

ground,  (fa)  Geometry  for  antenna  tilt  angle. 


Patterns  will  be  shown  for  the  antennas  tilted  through  sense 
angle  as  shown  in  fig.  29b.  note  that  the  antennas  arc  assumed 
to  be  tilted  directly  toward  tie  positive  x  axis.  When  an  antenna 
is  tilted  through  some  angle  «t  the  location  and  orientation  of  the 
infinitesimal  sources  representing  the  antenna  are  changed,  but  not 
their  magnitude  or  phase.  This  is  a  result  of  calculating  the  antenna 
currents  assuming  the  antenna  to  be  in  free  space.  Clearly  as  Oj-  in¬ 
creases  the  a/2  ground  plane  becomes  less  effective  in  shielding  the 
antenna  from  the  lossy  earth,  and  our  "free  space  assumption"  becomes 
less  valid. 


in  Figs.  30  -  33  the  far- field  power  patterns  for  antenna: 
are  shown.  Patterns  are  shown  for  6*  =  0C,  30°,  end  50=,  and 
the  xz  and  yz  planes.  The  patterns  are  normalized  so  that  th< 
into  each  antenna  is  one  watt,  and  so  that  the  peak  power  rad 
antenna  1  corresponds  to  0  db. 


in  both 
e  power 
iatea  oy 


At  low  elevation  angles  and  ot  =  0°  the  3-1/2  turn  MTL  {antenna  2 
is  several  dB  below  the  \f 4  monopole  over  four  */ 4  arms  (antenna  1). 
However,  when  fjt  =  60°  and  for  low  elevation  angles  near  the  positive 
x  axis  antenna  2  is  several  dB  superior  tG  antenna  1.  In  Figs.  32  and 
33  the  12  turn  picket  fence  MTL  is  seen  to  have  a  very  un symmetrical 
pattern  in  the  xz  plane.  This  is  largely  due  to  the  antenna  being 
fed  unsynrnetrically.  If  a  synuetric  pattern  ;s  desired,  a  symmetric 
feed  is  required.  Comparing  Figs.  32  and  33  we  see  that  the  4  inch 
plate  above  the  MTL  has  almost  no  effect  on  the  far-field  pattern. 
Antennas  3  and  4  exhibit  behavior  similar  to  antenna  2  in  that,  at 
low  elevation  angles  and  m*  =  0°  antenna  1  is  several  a3  superior, 
while  at  larger  the  12  turn  picket  fence  f-TL  may  be  superior. 

MTL  antennas  have  been  investigated  experimentally,  and  were  found  tG 
have  patterns  superior  to  the  quarter  wave  monopcle  over  fou’r  quarter 
wave  arms  at  low  elevation  angles  for  the  case  o*  =  60c.  These  measure 
merits  will  be  presented  in  a  future  report. 


It  is  possible  to  convert  the  relative  power  in  dB  values  rrom 
Figs.  30-33  directly  to  the  absolute  magnitude  of  the  field  strength 
in  volts/meter.  This  can  be  done  via  the  relationship 

p -?n 

(33)  !£j  =  (2.25  x  10  *  *  10  '  volts/meter. 

Here  P  is  a  relative  power  value  taken  from  Figs.  30-33  and  will  be 
less  than  or  equal  to  zero.  Equation  (38)  assumes  1  watt  ot  real 
power  input  to  the  antenna. 
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CURVES  FGR  VARIOUS  TILT  ANGLES  * 

Fig.  30.  Power  patterns  for  antenna  1  versus  tilt 
(a)  Pattern  in  the  xz  plane. 


Power  patterns  for  antenna  1  versus  tilt 
(b)  Pattern  in  the  yz  plane. 


CURVES  FOR  VARIOUS  TILT  ANGLES. 

Pig.  31.  Power  patterns  for  antenna 


f-ig.  31.  Power  Patterns  for  antenna  2  versus  til 
(b)  Pattern  in  the  yz  plane. 


a  nter 


Power  patterns  for  antenna  3 
(a)  Pattern  in  the  xz  plane. 


VARIOUS  TILT  flNGLr.S- 

33.  Power  patterns  for  antenna  3  versus  ti 
(b)  Pattern  in  the  y i  plane. 


VI. 


SUMMARY  AND  CONCLUSIONS 


In  this  report  we  have  had  two  basic  aims  directed  toward 
accomplishing  the  objective  of  this  investigation  as  stated  in 
Section  I.  These  aims  will  be  sunmarized  in  this  section.  The 
first  aim  was  to  obtain  a  qualitative  indication  as  to  how  the 
presence  of  a  lossy  earth  affects  the  far- field  VHF  power  pattern 
of  an  antenna  located  in  air  above  the  lossy  earth.  Our  second  aim 
was  to  describe  a  computationally  straightforward  metnod  for  determ¬ 
ining  the  far- field  pattern  of  an  antenna  located  in  air  near  a 
lossy  earth. 

To  obtain  a  qualitative  indication  of  the  lossy  earth's  effect 
it  was  decided  to  study  in  detail  the  infinitesimal  horizontal  and 
vertical  electric  current  elements  instead  of  specializing  to  a 
specific  antenna  geometry.  First,  computationally  useful  expressions 
for  the  fields  radiated  by  the  current  elements  were  presented  and 
it  was  shown  how  these  expressions  could  oe  interpreted  in  terms  of 
space  waves  and  surface  waves. 

Space  wave  fields  have  the':r  maximum  radiation  away  from  the 
earth's  surface,  are  characterized  by  a  1/R  radial  dependence  in  the 
far-field,  and  vanish  at  the  earth's  surface.  Surface  waves  have 
their  maximum  radiation  at  the  earth's  surface,  and  are  characterized 
by  a  1/R2  radial  dependence  in  the  far-field  for  frequencies  above 
approximately  100  MHz  and  typical  ground  constants.  Thus,  the  surface 
•wave  is  only  inportant  at  very  low  elevation  angles. 

Far-field  power  patterns  are  shown  for  horizontal  and  vertical 
elements  located  between  0  and  2.\  from  the  earth's  surface.  It  was 
found,  as  expected,  that  vertical  elements  are  better  radiators  at  low 
elevation  angles,  while  horizontal  elements  are  better  radiators  at 
large  elevation  angles.  Elevating  an  element  from  the  earth’s  surface 
was  found  to  generally  increase  the  radiation  at  low  elevation  angles. 

Far-field  power  patterns  were  also  shown  for  the  horizontal 
and  vertical  element  located  ,\/4  above  the  earth  and  tor  varying 
ground  parameters.  It  was  found  that  near  200  MHz  varying  the  ground 
conductivity  from  0.1  to  0.0001  mho/meter  has  little  effect  on  the 
element  patterns.  However,  varying  the  relative  permittivity  from 
2  to  32  increases  the  radiation  at  low  elevation  angles  for  vertical 
elements,  but  decreases  it  for  horizontal  elements. 

Because  in  the  far-field  the  space  ware  field  varies  as  1/R 
and  the  surface  wave  field  varies  as  1/R2  some  error  will  be  intro¬ 
duced  in  extrapolating  a  pattern  from  one  radial  distance  to  another 
if  tha  usual  1/R  field  dependence  is  assumed.  At  frequencies  greater 


than  160  MHz  and  typical  ground  constants  it  was  found  that  extrapo¬ 
lating  from  R  =  1  km  to  greater  radial  distances  can  be  done  with 
little  error  since  the  surface  wave  at  these  frequencies  attenuates 
rapidly  with  distance  from  the  source. 

Determining  the  far-field  pattern  of  an  antenna  can  be  thought 
of  as  consisting  of  two  steps.  First  the  current  on  the  antenna  must 
be  determined.  If  the  antenna  is  located  near  a  lossy  earth  the 
exact  solution  will  involve  an  integral  equation  which  contains 
the  Sontnerfeld  integral,  and  will  be  very  difficult  and  time  consuming 
to  solve.  Our  approach  was  to  assume  that  tne  antenna  current  was 
unchanged  by  tne  presence  of  the  lossy  eartn,  and  thus  it  could  be 
found  using  available  computer  programs  for  finding  antenna  currents 
in  free  space.  This  approximation  is  best  applicable  when  the  antenna 
is  shielded  from  the  lossy  earth  as  by  a  cavity  or  a  small  perfectly 
conducting  ground  plane,  or  if  the  antenna  is  physically  removed 
several  wavelengths  from  the  earth. 

Regardless  of  how  the  antenna  current  was  obtained  the  problem 
still  remains  to  find  the  radiation  pattern  of  a  now  known  current 
distribution  in  the  presence  of  the  lossy  earth.  The  exact  solution 
involves  an  integration  over  the  antenna.  The  integrand  will  contain 
the  Sommerfeld  integral  and  again  will  be  difficult  and  time  consuming 
to  evaluate.  Our  approach  was  to  replace  the  continuous  current 
distribution  on  the  antenna  by  a  finite  number  of  infinitesimal 
electric  current  elements.  The  fields  of  the  individual  elements 
can  be  evaluated  using  expressions  developed  by  Norton[10],  and  then 
surrened  to  yield  the  antenna  pattern.  The  accuracy  of  this  approach 
is  discussed  in  Appendix  II. 

The  above  approximations  were  applied  to  calculate  the  far-field 
power  patterns  of  several  antennas  located  near  a  lossy  earth. 
Specifically  we  compared  the  pattern  of  a  quarter  wave  monopole 
above  four  quarter  wave  arms  to  the  patterns  of  a  few  electrically 
small  multitum  loop  (KTL)  antennas  placed  over  a  half  wave  diameter 
perfectly  conducting  ground  plane. 

None  of  the  antennas  considered  had  an  overall  superior  far- 
field  pattern.  However,  for  large  tilt  angles  (e.g.,  ot  =  60°)  with 
the  antenna  tilting  generally  in  the  direction  of  the  observer,  the 
theoretical  calculations  tend  to  indicate  a  higher  level  of  radiation 
at  low  elevation  angles  from  the  MTL  than  from  the  quarter-wave  dipole 
(e.g..  Fig.  31a  and  Fig.  30a  respectively).  This  geometry  was  of 
particular  interest  in  this  investigation  since  it  depicts  the  situation 
wherein  an  MTL  mounted  on  a  projectile  would  impact  the  earth  at  a 
large  tilt  angle  and  transmit  back  to  the  area  from  which  it  was 
1 aunched. 


In  addition,  since  the  theoretical  calculations  described  in 
this  study  were  made,  the  patterns  of  several  MTL  configurations 
superior  to  those  in  Section  V  as  well  as  the  quarter  wave  monopole 
over  four  quarter  wave  radial  arms  have  been  investigated  experimentally. 
These  measurements  indicate  an  even  greater  superiority  of  the  MTL 
over  the  quarter  wave  monopole  for  the  geometry  of  particular  interest 
described  above.  This  is  most  likely  due  to  the  fact  that  when  et  2  60° 
the  radials  do  not  shield  the  monopole  from  the  earth  and  the  theoretical 
calculations  can  be  expected  to  be  in  error.  Apparently  this  error  tends 
to  cause  a  theoretical  prediction  in  field  intensity  greater  than  that 
which  is  observed  in  practice  for  the  monopole.  Since  in  a  communication 
system  the  important  factor  is  the  minimum  assured  power  level,  and  our 
measurements  showed  this  minimum  level  to  be  several  d3  better  for  the 
MTL  than  for  the  monopole,  we  conclude  that  the  MTL  has  better  radiation 
characteristics  for  the  application  of  primary  interest  in  this  work. 

An  area  of  future  work  would  be  to  determine  better  MTL  configurations. 
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APPENDIX  I 


In  this  appendix  two  computer  programs  will  be  described.  They 
are  DIPOLE  and  WERFC  and  form  the  basis  for  all  pattern  calculations 
macr  in  this  report. 

Figure  34  is  a  listing  of  subroutine  DIPOLE.  Included  is  a 
brief  description  of  the  program  and  definitions  for  the  12  input 
and  3  output  parameters.  Of  these  EXCIT,  £X,  EY,  and  EZ  are  complex, 
and  the  remainder  are  real.  The  program  is  an  extension  of  equations 
developed  by  Norton[10]  to  include  electric  current  elements  of 
arbitrary  position  with  respect  to  a  fixed  coordinate  system.  In 
Fig.  35  is  shown  an  arbitrarily  orientated  electric  current  element 
located  at  the  rectangular  coordinates  (XS,  YS,  ZS)  meters.  The  xy 
plane  is  assumed  to  separate  a  lossy  half-space,  described  by  a 
relative  dielectric  constant  £R  and  a  conductivity  SENU  expressed  in 
emu  (1  mho/meter  =  10“^  emu),  from  air  or  free  space.  The  lossy 
half-space  is  located  at  z  <  0.  The  source  is  assumed  to  have  an 
e~j“t  time  dependence  with  frequency  FKC  in  kilohertz.  It  is  desired 
to  evaluate  the  field  at  a  point  described  by  the  cylindrical  coordin¬ 
ates  (RF,  PF,  ZF)  with  RF  and  ZF  expressed  in  meters  and  PF  in 
radians.  The  orientation  of  the  current  element  is  described  by  THETA 
and  PHI  in  radians  as  shown  in  Fig.  36.  Here  the  x'y'z'  rectangular 
coordinate  system  is  the  xy-:  system  of  Fig.  35  translated  to  the 
source  point  (XS,  YS,  ZS).  The  subroutine  returns  EX,  EY,  and  EZ 
which  are  the  x,  y,  and  z  components  respectively  of  the  electric 
field  in  volts/meter.  NBUG  enters  the  subroutine  through  a 
CGMMON  statement  and  is  used  ir.  debugging.  If  NBUG  =  1  the  messages 
"ENTER  DIPOLE"  and  "EXIT  DIPOLE"  will  be  printed  when  the  subroutine 
is  entered  or  exited  respectively. 

DIPOLE  may  be  used  only  when  the  source  point  and  the  field 
point  are  in  the  air  half-space,  and  are  separated  by  at  least  a  few 
free  space  wavelengths.  The  program  has  one  difficulty.  Errors 
sometimes  occur  when  calculating  the  fields  very  near  (less  than  one 
degree)  the  vertical  axis.  This  difficulty  can  be  overcome  by 
interpolating  the  fields  at  either  side  of  the  vertical  axis  to 
obtain  the  fields  at  the  vertical  axis.  In  the  far-field,  the  fields 
are  varying  slowly  and  interpolation  may  be  done  with  high  accuracy. 

DIPOLE  uses  the  complex  function  subprogram  WERFC  in  evaluating 
the  ground  wave  attenuation  functions  of  Eqs.  (9c)  and  (15b).  Figure 
37  is  a  listing  of  WERFC.  Specifically 

_z2 

(39)  WERFC(z)  =  e  erfc(-jz),  z  =  x+jy. 
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SUBROUTINE  WEKFC  IS  USED. 

<t:«ses33fsssss:ts;ss5»;.ts595ssss*s-t«s5sss:;  »:99»»99;«39»t;3» 


Lt»i99C*9St««sr«St»59;SS9t5et«S9!'9  5S95tSSC9*SSS*a9*St;j59tS9S*tS959«9» 

COMPLEX? XI .Fx2.Fl .E2.X1 ,X2,X2,91 .V.G.EXCl i 

COHPLEXXK2,XJ,UtPP,HP,RVtF,Pl,H.  P.H  ,E2  V,  ER  V,  EP  V,  E2h  ,  ERH,  EPH,  Ex, 
2EY.E2 

COHPLEXIX,lY,I2tRE,EP,ERPH,EPPH,E2PH 

COHPLEXKERFC 

COMHONNBUC- 

IFINBUG.EO-l IWRJTE  I6tl  I 
1  FORKATI5X, 'ENTER  DIPOLE') 

XJ=CMPLXlO.O, 1.0) 

P 1=4. 0=A I AN (1.0) 

A=2S 

C=2.99BE3 

WAVE=C/IEKC=100n.0) 

XK  =  2  *0=P I /HAVE 


SET  UP  COORDINATES. 

XPF=RF*COSfPF I-XS 
YPF  =  RF  eSIN (PF  I-YS 
2Pf=2F 

RPF=SORT(XPF=XPF*YPF*YPF  ) 
PPF=AlAN2IYPF,XPr ) 

R=  SQK T <  RPFSRPF-. /PF  =2. "F  I 
XX=1.8ElH«Sf MU/FKC 
XK2=CSDRT IfK+XJ-XX UXK 
U=XK/XK2 


Fig.  34.  Computer  listing  of  subroutine  OIPOLE. 
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non  non  non  non 


R2*SC‘R  T  (P.PF*RPP,  (ZPF*  A  )*  (ZFF+A  J  ) 

R1-S0R1  (RPFaKPP*(/PF-AI*(ZPP-A>  ) 

SNC-ZPF/K 

SNCP--1  /PF«a)/P.2 

pp^aj*xk*r ?*u-u/; .0 

WP=PP*( l.Ct (A+/PF )/ <U«R2) )~~Z 

CSCP=SUkT  ( 1  .n-SiJO-S'JCP! 

rv=u*csort  ( i  .o-(ij-cscP)  =  (u-cscpn 

RV-1 SNCP-RV 1/ t$UCP*AV ) 

U=A1AN(  (EP+CSCP=CSCP)/XX) 

P  =  P  I  *K2 *CliS  l fS  i  /  I  XX -WAVS  ! 

PI-P-CCXP(XJ-B) 

W=4.0*?l/ i (1.0-RV )-(l .5-pv  I  j 
F  =  X  .0*XJ~CSCP.T  !  P  2  ~w  )  -wc  RFC (CSORT  (W  )  ! 

SNCPP=(Z?F-A)/R1 

CSCPP=S0R1 (1.0-.SNCPP*=2  J 

EXl=XJ*XKiRl 

EX2=XJ=:  <-R2 

E1=CEXP ic  xi 1/Rl 

E2=CEXP(£X2 1/R2 

Xl  =  ( 1-0/fcXl )-( 1 .  1/FX 1**2  ) 

X2=< 1.0/E X2 )-< 1.0/2 A2=*2> 

X3=CS0RT 1 1 .0-U**2*CSCP**2 1 

0E7ERH1NE  E2V 

EZV=CSCPP  =  *2*rl+P.V*CSCP  =  *2  =  E2 
E7V=EZV-»  1 1  .t-RV  }  s  ( 1  .0-t;*i:-Hi=*4*C SCP**2) *F  «£2 
6ZV=EZV-U-/.3=J2.0«SMCH-f  2/EX2 
E/V=EZV-ElrXl*n  .0-3.0*SNCr‘P) 

£ZV=£ZV-E2*X2* ( 1 .0-3.0'=SfJC=  ) 

£ZV=XJ*XK=EZV 
A1=CABS (fc  ZV*R/ ( 2 .0*XK  )  ) 

DETERMINE  EP-V 

ERV=SNCPP*CSCPP*rl+RV*SNtP-C$CP*£2 

ERV=ERV-CSCP-(l.O-RV)*U*X3*Fse2*U.o-0.S<U*»2*X3e*2+0.S*SF.,CP«*2-O. 

23/EX21 

ERV=ERV+5NCP*CSCP*  ( 1 .0-RV 1*E 2/EX2 

rRV=£RV-3.0*S.NCt,?*CS:.PP;;Xl  =E i+CSCP*U*X3* Cl. 0-RV) *2  2/ (2.0-EX2 1-3.0* 
2SNCP*CSf P5X2-EZ 

£RV=~XJ*XK*ERV 


DETE8H1NE  EPV. 

EPV=CHPLXCC. 0,0.0  ) 

RH=CX3-U*SNCP)/  C*3-*USJKCP  1 
BP=ATAN(CtR-CSCP**2»//x  1 
0=P 1 *XX*R2/ (COS !SP )=WAVc 1 
01=-0*CEX?I-XJ*PF) 

V=4. 0*01 /C1.C»RH 1**2 

G= 1 .0+Xj*CSG£T I P I *V 1 =k£RFC 1 CSOST (VI)  "  '  ‘  '  -  ~ 

DETERMINE  EZH. 

EZH-  SNCPP=CSCPR*F 1-RV*SNCP*CSCP*£2+CSCP*C1. 0-R V) *U3X3*F*EZ*C 1.0- 
2UsU*0.5«X3®*7-*o.$*SNCPaS2-C.5/Ex2 1 
eZH=FZH-SNCP-C\!  P-  { 1  ..‘,-RV )SE 2 /EX 2-3.0* SNCPPsCSCPP*Xl*£l 
fcZH=fcZH-CSCP*  ( 1  .O-RV  ;~’J-x3~E2/i  2.0*EX21  ♦3.0-5nCP*C5CP*X2=F2 
EZPH-XJsXR-tZH 
EZhsEZPH-COSIPPF 1 


Fig.  34. 
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ooo  ooo  ooviooo 


DETERMINE  £RH. 

ERH=SNCPP==2<£l-RV=SNCP=*2=E2+< i.O+RH )-G=E2-CSCP*=>2*U*«2-( 1 .O-Ry ) 
2'f*E2 

ERH=ERH-X1~ I 1.0-3 .0-CSCPP“s2  >-S 1  +  X 2- ( 1. 0-3. 0=5CSCP«s2 ) « ( 1 .0-U«2= 
2(1.0+RV)-U**2*il.3-<V  i~r )-£2 

ERH=ERH+U*»2SCSCP- -2- I 1 .0-RV )={ 1. 0- 1 . 0/6X2 )~ 1 F= { U==2=X3-SNCP=*=2 
3-1.0/EX2)  +  l..'/£X2)==2 
£RPH=— XJ-XK1 cKri 
ERH=ERPri=COS(PP= ) 


06 Tt RHINE  EPri. 


EPH=  f l-RH*  F2 ♦ ( 1  .O'Kh  )=•»{•?-(  i.  ')-i.O/£X  l)«l  1  /  EX  1 

EPH=EPrt+  u  ,(J-1  «0/t  X2  )°  I  1  .  (’-•>'  *2  -  (  J  .  f-.RV  l-U*  '  /■»  i  l.O-RVI  ->E  )  *f  //t  X ? 
EPH=f  PH+O.b -!)-»?=  I  1  .U-iV  )  =  (r  -  <i;=*2«-X3-  i.\CP*«2-  !.0/EX2)  +  l  .O/1  X2)  * 
2E2/EX? 

EPPH  =  XJ<-XK=EPH 
EPH=£PPH»S1MPPF  ) 

PLACE  PROPER  MAGNITUDE  AND  PHASE  ON  ALL  COMPONENTS. 


CONVERT  TO  RECTANGULAR  COORDINATES. 

I X=6XC I TsSJNCTHFTA )iCOS(PHl ) 
IY=EXCIT»SINITH3TA JSS1N (PHI ) 

I2  =  EXC  J  T<-CO$  (  THE  7  A  ) 
RE=Ix^erh+jy5;6p.phssimppf  )+I2=ERv 
EP=  1  X=£PH-1  Y-"*E PPn*COS  (?PFl 
6  2=  IX=E2K4|Ys£2PH-SiN(-’PF  HI 2 =E2V 
EX=RE*COS(PPr  )-rp-SI.\(-P=) 

£Y=RccS IN  I PPF I +f P=COS (SPF ) 

If  INB'JG.EQ.I IwKITE (6,21 
2  FORMATC  20X, *6X17  OIPULc'l 
RETURN 
END 


Fig.  34. 


57 


.  35.  Location  of  the  source  in  subroutine  UIPOlE 


a .  36.  Orientation  of  the  source  in  subroutine  D 


o  o  o  o 


CnMPiexRINCl IHNWERFC ( l ) 


WFRFCIZ!  eviL'iATES  KIZ)  =  EXP(-Z=«2)eERFC(-J=2)  USING  N.B.S. 
EOS.  7.1. ?3  A*‘D  7.1.29. 


CfVP!  FXKEP-FC.7,XJ,7Z,rPFCtA,B,C 
CnUHUNNBUG 

C p-Sh ( x  i  -  I f  XP ;  V  } XP  { -x  i  # /? .G 
S 1  NH  l  X  J  r  (  E  XP  I  x  >  -P  XP  f  -X  5  j  /  2.  0 

F  CXtY,N)=?.f.sx-2.0*X-C0S!!(.N«Y  >=COS  I  2. 0  =  X=Y  1 4N«S  INHI its  Y  }  =  S  IN  (2 .0*X  = 

2Y) 

G(X,Y,N  )=?.iRXsr.i><:h  (?;JY  )£<; in !  2.  0*X  =Y1  +N«S  1NHI  K-  y ) 5 COS (  2.0=i=Y ) 

IF  (NfiuG.FQ.l  imp:;,  jj 
1  FORMAT  I  SX  <  “-NT  t-P  RBFC'l 
PI =4.0‘ATAM( l.OI 
X J=CMPLXfO.O, 1 .0) 

Z7 =-XJ=Z 
X=S£ALI77) 

Y=  A  l  M  A  G I  Z  7  ) 

C  wen  E  <6,90)X,V 

IF  ( ABS  (X  1  .GT  , 3  ,Q.;iP.  AE»S(  Y  I.GT. 3.0) G0TO200 
T  =1  .0/1  l.C+.P7«!!-'X  I 

Fp  Fx  =1  .0-(.2S*H?c*.n?s:T-.2flAA9A7  3^T  =  e2+l.<.2JAl  374!*T**3-1  .AS315202 

27±Ts»A+1.06IAG?A?q5Tci5)'EXP I-X=X 1 
C=EfcFX-*(6XP(-X=X)/  I  ?.0SP  l  £x  )  1*1  I l.O-C0S(?.0=X=YJ  ) +  XJS=S  INI  2  .0*XSY  >  1 
AiCMPLX (0.0.0. 01 
1  =  0 

oniooN=  1,100 

B=  IEXPI-.25«iSNl/{fKN+A.0*X=X  n*IF(X,Y,N)  +  XJ»G(X,Y,Nl  ) 

A  =A+B 

C  WB|TE(6,90JA,B 

90  FrRMATI'.X.AtlA.B) 

IF  (N.LE.2  1GOKMGO 
Rl=ABSIkrAL ( B 1  ) 

R?=ASS (REAL (i) ) 

Sl=ABSl A1MAG!» )  ) 

S?=ABS(AIKAGI A) ) 

I F (Rl/ 1 P.l  +R2  I .LT .3.OE-6.AN0.  S1/ISI+S2)  .LT. 3.0E-6) G0T01 10 
1  =  0 

GOTOIOO 
110  1=1+1 

1  FI  I .GE .3 ) SOTO 120 
ICO  CONT INUE 

WBITE(G,!?0J 

130  FORMAT  I//SX,'  =  P.FC  FAILED  TO  CONVERGES//! 

120  A=2.0=EX!*I-Xix  K  A/P  I 
E°FC  =  l  .0- 1  A*r.i 
WFRFC=ERFCiCEYPI-7=Z) 


Fig.  37.  Conputer  listing  of  function  subroutine  WERFC. 
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GOT  0300 

200  erfc=cmplx<i .0,0.0) 

I  =0 

N«AX=CAfi$<ZZ> 

NMM=i 

oni50N=i  ,nxax 
A=EREC 

E  RFC  =E  R  r C  +.NUM  /  i  2 . 0-  22  =2  2  )  N 
IE  tfi.LE.l  )GOTfU50 
xa=cabsia: 

aE  =CABS ( £  RFC  ) 

RAT=ABS iXA-XE  ?/X E 
IE  (P.AT.LT.l  .0F-5  IGOTCifcO 
I  =0 

GOT 01 50 
160  1=1  +  1 

IFII.GE.2  1G0TO 170 
150  CONTINUE 

IE ( RA .  ,T.-.0F-3)HRIT£(6,1B0)RAT 
1B0  FORMA  a.'HFRPC  FAILED  TC  CONVERGE  RATIO  =  ’.E13.4/I 

170  WFRFf  •  <•  (  27  <■  SORT  (  P I  )  ) 

300  CONTI 

IF (NBUL  .1 )<JRITE (6, 190)2, WERFC 
190  FORMA M „a, *2  =  ',2E15.5,5X,'KERECIZ)  =  *,2E15.6) 

IF (NBUG.EO.l )UR IT  t (6 • 2) 

2  FORMAT120X, 'EXIT  WERFC* I 
RETURN 
ENO 


Fig.  37. 


60 


The  following  formulae  are  used[ll] 


(40)  &  z  ez  erfc(z)  *  1  +  \  (-l)m  i^3-  '!2^.1-) 

...  ^  o  Hi 


m=l 


(2z2) 


for  ( x |  >  3.9  or  jyj  >  3.0,  and 


-x 


(41)  erf (x+iy)  =  erf(x)  +  [(l-cos(2xy))  -s-  jsin(2xy)] 


o  2  -n^/4 

+  |e"x  l  —p — o  Cfn(x»y)  +  J9n(x,y) J 
n=l  n^+4x^  n  n 


otherwise.  In  £q.  (41) 

(42a)  erfc(z)  =  erfc(x+iy)  =  l-erf(z) 

(42b)  fn(x,y)  =  2x-2xcosh(ny)cos(2xy)+nsinh(ny)sin(2xy) 

(42c)  gn(x,y)  =  2xcosh(ny)sin(2xy)+nsinh(ny)cos(2xy). 

NBUG  enters  this  function  subprogram  through  a  COMMON  statement 
and  is  used  in  -tebugging.  If  NBUG  =  1  the  messages  "ENTER  WERFC” 
and  "EXIT  WERFC"  will  be  printed  when  the  function  subprogram  is 
entered  or  exited  respectively.  Z  and  WERFC(Z)  will  also  be  printed 
if  NBUG=1 . 
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APPENDIX  II 


In  this  appendix  far  field-patterns,  evaluated  using  subroutine 
DIPOLE  described  in  Appendix  I  and  approximations  described  in 
Sections  II  and  IV,  will  be  compared  with  patterns  from  the  literature. 
In  Figs.  38-40  we  compare  far-field  patterns  of  infinitesimal  electric 
current  elements  calculated  using  subroutine  DIPOLE  with  patterns 
presented  by  Jordon  and  Balmain[14].  The  ground  parameters  and 
frequency  for  these  patterns  are  identified  by  the  parameter  n  and 
by  as  follows: 

(43)  n  =  x/Ej. 
where 

18  x  10^  o,.  . 

(44)  v  =  _ mho/met-.r 

w  x  f  MHz 


Each  figure  shows  curves  for  n=l  and  n=100,  and  the  curves  are 
normalized  so  that  the  n=100  curves  have  maximum  field  strength  of 
1.0.  Figures  38-40  are  for  a  vertical  element  at  the  earth's  surface, 
a  vertical  element  a/2  above  the  earth,  and  a  horizontal  element  >./2 
above  the  earth  with  the  pattern  in  the  yz  plane  (see  Fig.  4) 
respectively.  In  addition  our  calculations  have  been  found  to  have 
the  proper  rotation  symmetries. 

Next  we  will  investigate  an  extended  source,  namely  a  vertical 
half-wave  dipole  with  its  center  height  one  wavelength  above  a 
lossy  earth.  The  frequency  is  20  MHz  and  field  patterns  are  shown 
for  the  elevation  plane  (see  Fig.  4).  Assuming  the  antenna  to  be  in 
free  space  we  approximate  the  antenna  current  as  being  sinusoidal. 
Using  the  trapezoidal  rule  and  the  methods  of  Section  IV  the  antenna 
was  represented  by  11  vertical  infinitesimal  electric  current 
elements.  The  fields  of  these  11  Infinitesimal  dipoles  were 
determined  using  subroutine  DIPOLE  and  then  summed  to  yield  the 
patterns  shown  in  Fig.  41.  Patterns  calculated  by  Hom[15]  are  also 
shown.  Fig.  41a  is  for  c  =  0.1  mho/meter  and  cr  =  30,  and  Fig. 

41b  is  for  o  -  0.01  mho/meter  and  er  =  10.  Our  patterns  and  those  cf 
Horn  have  been  normalized  so  that  the  peak  field  strength  of  the 
patterns  shown  in  Fig.  41a  are  unity.  Note  the  excellent  agreement  cf 
the  location  and  magnitude  of  maximums  and  nulls  in  Fig.  41. 
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JORDON  a  BALMAIN 
PRESENT  RESULTS 


JORDON  a  BALMAIN 
PRESENT  RESULTS 


ffi 


HORN 


center  height  ,\  from  a  lossy 

(a)  o=0.1  mho/meter,  e  =30. 

(b)  o=0. 01  mho/meter,  e  =10. 
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